Diffusion of a Deformable Body in a Randomly Stirred 

Host Fluid 

Moshe Schwartz, Gady Frenkel 

Raymond and Beverly Sackler Faculty of Exact Sciences 
School of Physics and Astronomy 
Tel Aviv University, Ramat Aviv, 69978, Israel 

We consider a deformable body immersed in an incompressible liquid 
that is randomly stirred. Sticking to physical situations in which the body 
departs only slightly from its spherical shape, we calculate the mean squared 
displacement and the diffusion constant of the body. We give explicitly the 
dependence of the diffusion constant on the velocity correlations in the liquid 
and on the size of the body. We emphasize the particular case in which the 
random velocity field follows from thermal agitation. 

I. INTRODUCTION 

Systems of deformable objects immersed in a liquid are very common in every day life. 
Milk and blood, for example, are such composite systems. Milk can be viewed as an emulsion 
of fat globules in water while blood is a suspension of cells ( that have some rigidity) in water. 
The physical description of the set of objects present in a given liquid involves the location 
of the objects, their shape and in some cases the strains on the objects or any other fields 
that are needed to describe the objects in addition to their location and shape. The actual 
solution of such systems is extremely difficult because each object interacts with itself and 
with the other objects via hydrodynamic interactions. Hence, we are facing a many body 
problem with the additional complication, that each object is not described by a single 
degree of freedom ( its center of mass) but actually by an infinite number of degrees of 
freedom, where all the deformation degrees of freedom, corresponding to all the objects, 

1 



interact. The situation is simplified a little if the deviation of the objects from spherical 
shape remains small This happens when the agitation of the host liquid is not too 

strong and when the density of the objects is not too high ( Close packing would cause finite 
deformations, although perhaps still treatable within the small deformation approximation). 
Our final goal is to obtain the response of the composite system to a given velocity field 
imposed on the liquid. The velocity field we have in mind may be fixed in time like simple 
shear or randomly fluctuating in time and space. Even in the first case the velocity field 
experienced by each object separately must have a random part due to the random passage 
of other objects near by. The present article is the first in a series that deals with this general 
problem and it concentrates, within the small deformation approximation, on the diffusion 
of the center of mass of a deformable body in the presence of a random external velocity 
field, imposed on the liquid. It is possible to consider the center of mass separately from the 
deformation degrees of freedom, that will be discussed in a future publication, because as 
will be shown in the following, they are decoupled in the small deformation approximation. 
We study the mean square displacement (MSD) of an object as a function of time for a 
general random velocity field that has given correlations in space and time. For long periods 
of time the MSD is usually linear in time, enabling us to discuss it in terms of a diffusion 
constant that depends on the size of the object, R, and the correlations present in the liquid. 
There are, of course, cases where the MSD does not behave linearly at long times and the 
method we develop here is quite capable of dealing with those cases too. Our main concern 
is in objects in which the state of lowest energy is of spherical shape. This is the case for 
deformable objects dominated by surface tension [Q. Our results results concerning the 
motion of the center of the object will hold also for cases where the shape of lowest energy is 
nearly spherical (for example a body with bending energy and spontaneous curvature 
close to that of a sphere with the same volume). 

The paper is organized in the following way: In section ^ we define the system under 
consideration, discuss its properties, and construct the basic equations. In section ^ we 
derive the equation of motion for the center and the MSD in integral and differential forms. 
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In section |^ we discuss the results and demonstrate their use for different kinds of noise 
reahzations. The specific case of thermal agitation is considered in section |V]. In appendix 
^ we construct the velocity correlations for the case of thermal agitation and in appendix 
we consider the physical conditions under which the small deformation approximation is 
valid. 



II. THE SYSTEM 

Consider a single deformable body immersed in a host fluid. The system is chosen to 
have the following characteristics: 



1. The host fluid is incompressible, V ■ v = 0. In addition, we assume that the Reynolds 
number is small so that the Stokes approximation is applicable. 

2. The body is characterized by an energy that depends on its shape. The shape of 
minimum energy is a sphere. The shape of the surface of the body is described by the 
equation: ilj{r) = where ipl-f) is a scalar three dimensional field (Fig. [l|).( Although 
our derivation considers only objects for which the shape of minimum energy is a 
sphere, all the conclusions concerning the MSD carry over to cases where the shape of 
minimum energy is nearly spherical). 

3. A deformation of the body induces a force density on the host fluid. As a result, the 
fluid's velocity is given by the sum of Vext, which is caused by external sources, and 
v^, which is induced by the body. 

4. The external velocity, Vext, is random and is chosen to have zero average and known 
correlations. It is convenient to define the external velocity in terms of its spatial 
Fourier transform as 
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where u{q) is a general vector field and the subscripts denote Cartesian components. 
This definition implies only that the fiuid is incompressible and in any other way is 
general. Next, invariance under translations in space and time and under rotations 
yields the form of the correlations of the velocity, 

{ui{q, t)) = and 



where 6im is the Kronecker delta and 6{) is the Dirac delta function and is a general 
function of q and |t2 — If the random velocity field is characterized by a length scale 
and a characteristic time scale then it is convenient to write it as 0(^9, |t2 ~ til/'^)^ 
where ^ and r are, respectively, the correlation length and the memory time scale of 
the external velocity . 

5. The surface elements of the body are carried by the host fiuid 0, i.e. each surface 
point moves according to 



6. We assume that the external velocity is weak enough to cause only minor shape fiuc- 
tuations of the body. 

We will be interested in the following in the Mean Squared Displacement (MSD) of the 
center. Since the body is deformable the definition of its center is not unique. For periods of 
time shorter than r the result depends on the definition of the center. It turns out, however, 
that the value of the MSD at longer times does not depend on the specific choice, because 
for long times the MSD (according to any reasonable definition) is much larger than the size 
of the body. Therefore, the results for the diffusion constant are general and do not depend 



{ui{q,ti)um{p,t2)) = 6imS{q + p)(j){q, 1^2 - ti\), 



(2) 



(3) 
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on the specific definition of tlie center wliicli will be determined later. In cases where the long 
time dependence of the MSD is not linear, it is still tending to infinity with time, so that 
again the specific definition of the center does not matter. 



Following the line of derivation of Schwartz and Edwards equation may be 

turned into a continuity equation for ip, 

^ + {ve.t + v^) ■V^^J = 0. (4) 

Consider a deformable body, carried by the host fluid in such a way that at any instant it is 
nearly spherical. Its state can thus be characterized by the position of its center, ro(t), and 
a deformation function that describes the shape by the equation 

V;(f,t)^| + /(fi,t)-l = 0, (5) 

where p = \ f—fo\ is the distance of the surface from the center in the direction of the solid 
angle fl and R is the radius of the body when not deformed. The deformation function / 
can be expanded in spherical harmonics, = J2'iZoY.m=-i fi,mi't)Yi,mi^)- The center 

of the shape, ro(t), is defined as that point around which fimif) = 0. 
Equations (H) and (|^) lead to a linear equation for each fi ^, 

9fi 1 

—g^ + >^lflm + ^[p-{Vext-ro)]lm = 0, (6) 

where p is a unit vector directed outwards from the center in the direction of Q, and 

[p ■ (vext - ^o)]im = j dn[p- [t;e.t(ro + R{1 - f)p) - Fo] Y*^{n)} . (7) 

The eigenvalues A/'s characterize the decay of a slightly deformed sphere into a sphere in the 
absence of the external velocity. Different physical systems are characterized by different 
sets of A;'s. However, it is obvious that the decay must depend only on / because of the 
spherical symmetry. Examples of systems for which different sets of A; have been calculated 
include: a droplet with a surface tension and equal viscosities inside and outside ^ and a 
droplet with surface tension for a viscosity much higher inside the droplet than outside . 



Other systems for which the following results are applicable to, in the small deformations 
approximation, include a droplet with a bending energy 0, a droplet with a bending energy 
and in-plane dissipation ||T0| and a droplet with both surface tension and bending energy 
Our following discussion is therefore general and not limited to one specific system. 



III. DERIVATION OF THE MSD 

Equation implies that in order that fi^m stays zero for all times we must have as an 
equation determining the location of the center 

[P ■ {Vext - = , m = -1, 0, 1. (8) 

For Z 7^ 1 it is clear that tq can be dropped from the last term on the left hand side of 
Eq. (^. Therefore f{Q,t) is linear in Vext (for long enough times the initial deformations 
have already decayed). Consequently we can always drop, for small enough Ve^t, f in the 
argument of Vext on the right hand side of Eq. (|^) (The physical conditions for which 
this approximation is valid are discussed in appendix |Bp. This results in decoupling of the 
deformation degrees of freedom from that of the center of the sphere. The equation for the 
motion of the center can thus be given, using linear combinations of Yi^^, in vector form as 

J dn p {p-f^ = j dVt p {p- Vext{ro + Rp)) . (9) 

We integrate the left hand side of the above and express the external velocity in terms of 
its Fourier transform on the right hand side to obtain 

^fo = [dn [d'qpip- vUq, t)) e-^^^-(''»+«^). (10) 
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We use the partial waves expansion [|T2|,ITE 



-iq-{Rp) 



CO / 



E {-i)h7^jMR)YL{^,)YU^), (11) 

1=0 m=-l 

where VL and VLq are the solid angles in the directions of p and q respectively and ji is the 
spherical Bessel function of order /. We integrate over f2 and obtain 

'ro = ^j dqe-"^-'^'^ (^Jo(gi?) +j2{qR)A^ Vext{q,t). (12) 



The matrix A{q) is given by 



A, = —S., + U, - ^ ) . (13) 



It may seem that A on the right hand side of Eq. (p!2[) mixes directions. However, the 
bracketed term in Eq. (|l^) is just a projection operator on the transverse direction. The 
external velocity is incompressible and hence already transverse. Consequently, this term 
acts as a unity operator, 6ij, and Eq. (|1^) leads to 

r-o = J dqe-"^-'^° {UqR) + j2{qR)) v^M, t). (14) 



Equation (14) is the explicit equation of motion for the center of the body. In the limit 
R ^ the approximation, tq = Vextifoii) is obtained. Note that this equation is general 
and describes the motion of the center for any given (small enough) external velocity field. 

Next, we calculate the MSD, ((Afo) ), clS cl function of the elapsed time, t. Consider a 
specific realization of the external velocity field. 

The Displacement of the center is given by the trivial equation 

Afo(t)= f'r^{t')dt'. (15) 

Hence, the MSD is given by 

((Afo(t))') = f dh f dh (fo(ti) ■ Fo(t2)) . (16) 



The correlations of the external velocity given in q space are obtained from Eq.(|l]) and @) 
. Assuming the decomposition [111,0 



;t^e.t.(gi,tl)t^e.i,(g2,t2)e-^«"^l-^^»(*l)e-''-^«(*^) ) = 

{vextMM) Ve.t,{q2M)) (e'^'T-^^^'^^ e''''^'-^°^''^) (17) 
and in addition that the distribution of Aro(^) is Gaussian, i.e. 

-ig-Aro{t)\ ^ g-4({Aro(t))2)^ ^^g^ 



we obtain 



{3o{qR) + J2M)'E (l - ^) 1^2 - til). (19) 

The only term that depends on angel is 1 — q^/q^- Performing the angular integration 
/ dVlq (1 — q]^/q^) = 8tt/3, and summing up the three terms we obtain, denoting the MSD 
by Fit), 

F{t) = 16n / dt' / g^^ge- t') (UqR) + j2{qR)f [t - t'). (20) 
JO Jo 



Eq. (^) can be turned also to a differential equation. Differentiating Eq. ( pOD twice we 
obtain 



/■OO 2 

F(t) = IGtt / q^dqe--^'-'U{q, t) {jo{qR) + 32{qR)f ■ (21) 
The initial conditions are 

F(0) = (22) 
and 

F(0) = 0. (23) 

The latter condition is valid in cases where the correlation function, (j){q, t), is finite at t = 0. 
The only exception is the case of white noise, where one must carefully check the result of 
the first differentiation and determine F{0). (Actually (^) is always correct, because any 
noise that is of physical origin must be correlated in time. The widely used white noise is 
just a very useful idealization of the real situation, that will result in -F(O) = and F{S) 
having a value that is not small for rather small 5's). The advantage of the differential form 
is that its numerical solution can be easily obtained by advancing F{t) in time. Note that 
equation (^) above is not restricted to cases that can be described in terms of a diffusion 
constant. 
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IV. PROPERTIES OF THE MSD 



The random velocity field may be caused by thermal agitation which is an equilibrium 
phenomenon or by non-equilibrium process such as mechanical stirring. While Eq. (pTD can 
supply, by numerical solution, the MSD for any velocity correlation, there are families of 
velocity correlations in which at least part of the solution of (^) or ( PU[ ) can be obtained 
analytically, rendering the process of solving for the MSD much easier. The simplest case is 
where the correlations are white in time, namely (piqjt) = (f){q)6{t). In those cases F{t) is 
linear at all times, F{t) = 3Dt, where D is the diffusion constant and Eq. that is an 
equation for the function F{t) is replaced by an explicit expression for the diffusion constant 

Q-jr roo 

D = - qUq miMqR) +j2{qR)f. (24) 
6 Jo 

A family of correlations that is a simple extension of the above, where it is quite easy to 
see what is happening, is defined by (f>{q,t) = (j){q)G{^), were G is a function that decays 
when its argument becomes of order 1. It is clear from Eq. ( p3D that for short times the 
MSD must behave as while for long times it must be linear in t, since for t » t the 
time dependence cannot be distinguished from white noise. The function 0(g), will naturally 
have a cut-off factor g{qC,), where ^ is the correlation length. Clearly, the correlation length 
cannot be expected to be smaller than the distance between the particles of which the fluid 
is composed and not larger than the size of the system. 

The MSD depends, of course, on the ratio 7 = j. Generally speaking, as 7 increases 
the slope of the MSD and particularly the diffusion constant decreases. This is due to the 
fact that as 7 increases, different regions of the surface become less correlated and move in 
different directions. In the limit 7 ^ 00, the movement of the center ceases and F{t) is 



always zero. In the limit 7 — 0, the bracketed Bessel term in Eq. (^) and (21) can be 
replaced by unity (since lim^/j^o (joiqR) + 32{,qR)) = 1). A close inspection of the derivation 
reveals that this limit produces the same MSD equation as the equation for the approxima- 
tion ro = Vexti'f^o)- I-e. the latter approximation is accurate for an infinite correlation length, 
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or point particles. 



In the following we will consider the dependence of the diffusion constant on the size of 
the object R. Consider a correlation such as (p{q,t) = C 6{t) {q^)°'g{qC,), where g is a. cutoff 
function and g{0) > 0. Note that as discussed above the results that will be obtained here 
for the diffusion constant hold true also for a finite correlation time. We insert the above 
correlation function into equation (|20|), then substitute qR with u, and obtain 



^ = 3R^ I "'^"^(|")boM + J2M]^. (25) 
In the limit R/^ —>■ oo we distinguish between two cases: a < 1 and a > 1. Since the large 
u dependence of jo{u) + j2{u) is proportional to cos{u)/v? we find that 
-^r^ for a < 1 

D oc I ^ (26) 
^ for a > 1. 

In the opposite limit R/^ — > we find that regardless of a 
C 

e 

(Note here that we have written the power law dependence of 0(g) as but having other 
dimensional constants in the model may make C depend on ^, so that Eqs. ( pG] ) and ( p7[ ) 
may be considered only as equations that yield the dependence of D on the radius R). 



D oc -. (27) 



The dependence of the diffusion constant on the decay time scale can be also deduced 
when the correlation function is separable (i.e. the second family). Using simple dimen- 
sional analysis, Eq. (|20|) leads to the conclusion that the diffusion constant is linear in r ( 
in addition to the possible dependence of C on r): D ~ Ct. 



There is another class of velocity correlations that is not separable but allows the calcu- 
lation of the long time behavior of the MSD. This class is defined by a scaling form of the 
velocity correlations, 

(P{q,t) = Cq~y {Tqt") , (28) 
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where T and C are dimensional constants, / is a function with a finite decay length and 
a < 3. The solution is obtained by assuming that F{t) = Af^. The integrand in Eq. ( |2TD 
has in it two functions, each cutting the integral off at different value of q that is a function 
of t. The dominant cut-off at large times t, is the one that cuts the integrand off at smaller 
g's. What remains is just a scaling argument that leads from Eq. (^) to the following 
result: 



(29) 




where 

4/(5 -a) ifz>>2/? 



(30) 

2 + (a - 3)/5 if z> < 2f3 



Note that Eq. (29) results from the fact that even if F{t) = for large t the leading behavior 



of F{t) is still linear. Note also that in Eq. (30) the two options have to be evaluated first in 



order to check which of the conditions applies. An explicit equation for the prefactor A can 
also be easily obtained. We solved Eq. (|2T|) numerically with a = 5/3 and /? = 3/5. The 
long time dependence of F(t) is depicted in Fig. (^. We see that the long time dependence 
is given by F(t) oc t^'^, which is exactly the result predicted by Eq. (|30D . (Note that in this 
case u = 2/5). 



V. THERMAL AGITATION 

Of particular interest is the case where the fiuctuations in the velocity field are due to 
thermal agitation. We describe the effect of temperature by a scalar potential ip and a vector 
potential A, that fluctuate, have zero average and local correlations in space and time. Both 
give rise to a force density fleld 

F = -Vip + V X A, (31) 
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that generates in its turn the fluctuating velocity fleld in the hquid. Since the velocity fleld 
is divergence-less, the scalar potential affects only pressure. Hence, the Fourier transform of 
the velocity fleld is given within the Stokes approximation by 

iq X Aiq, t) , , 

vm = - r^, 32 

where r] is the viscosity of the liquid. Since the correlations of A{f, t) are local in space and 
time, it follows that (j){q), deflned by Eq. @, is given by 

k<l) = ^, (33) 

where c is a dimensional constant. Dimensional analysis reveals that c must be proportional 
to KbT/t] (with a dimensionless proportionality constant). A detailed calculation yields a 
proportionality constant equal to (27r)^'^ (see appendix 0). The flnal conclusion is that for 
R larger than the inter-particle distance in the liquid |I6| , 



D = (34) 

Note that this result, for a liquid membrane that has liquid inside as well as outside, is 
different from the stokes result for a hard sphere. (We may expect Eq. ( P3| ) to hold only for 
g < 1/,^ where ^ is the inter-particle distance in the liquid but since R is expected to be very 
large compared to we are always in the situation described by Eq. (^61) with a = —2). 



The latter result is similar to the result for a polymer subjected to thermal fluctuations [17 
(with a different pre-factor). 

APPENDIX A: VELOCITY CORRELATIONS FOR THERMAL AGITATION 

We wish to determine the exact form of the velocity correlation function for the case 
of thermal agitation. Consider a system in which the random velocity fleld results from 
thermal agitation. The transversal part of the linearized Navier-Stokes equation reads: 

^ = -^I'vl + — , (Al) 
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where Vq is the Fourier Transform (FT) of the velocity field, v is the kinematic viscosity, F 
is the FT of the force density in the liquid, p is the number density of the particles and m 
is the mass of a liquid particle. 

We solve this equation under the condition that the force density is white noise in time and 
obtain, 

v],{tiH{h)) = e-^'^'\'^-'^\ {vlvl) (A2) 

where the last average on the right hand side is an equal time average. It is clear from the 
above that the velocities at different times are correlated, as opposed to white noise. It is 
possible however to consider effective white noise correlations by integrating the right hand 
side of (|A2|) over time and replacing then the decay function exp{—vq^\t2 — ^i|) by some 
a{q)6{t2 — ti), that will produce the same integral. This yields for the effective white noise 
velocity field, 

{u,{q,h)u,{p,h)) = ^5{t2-h){v\vi)^^, (A3) 

where we denote the effective velocity field by u and the real one by v. (v'^^v^^^ must be 
proportional to 5{'p + q) because of invariance to translations and to [5ij — because 



of incompressibility. Comparing with Eq. (|33D , we conclude that there is no additional 

dependence on q. 

Therefore, 

(«) = c5(t2-ti)[5u-^]. (A4) 

Now, we wish to relate the velocity to temperature. Considering that our continuous liquid 
is actually made up of discrete particles each having a mass m, we know that the total 
kinematic energy of the liquid is 3N^^f^ and as a result we find 

V«) = (A5) 
Expressing v in terms of its FT and integrating while keeping in mind that the number of 
degrees of freedom should be conserved and equal to 3A^ yields 

c' = — (A6) 
2mp(27r)3 ^ ^ 
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from which we can see that 

and the full correlation function reads 



{u,{lh)u,{pM)) = + ^^(^2 - h)[6,j - ^]^. (A8) 



APPENDIX B: VALIDITY OF THE SMALL DEFORMATION APPROXIMATION 

The validity of equation (pOD is limited by the approximation of replacing Vext{rQ + R{l + 
f)p) by Vextiff) + Rp) that has been discussed in section |ITI[ The approximation implies 
that we can, in the limit of small deformations, replace the velocity at the surface with the 
velocity on the undeformed sphere. To check the approximation, we expand equation (H) to 
the first nontrivial order in /, 

j dnYl^{n)p ■ [vextiro + Rp) + Rfi^, t) {p ■ V) Vextiro + Rp) - ^o] = 0. (Bl) 

The approximation is justified if the first order term is negligible with respect to the zeroth 
order term. A careful inspection reveals that this condition holds if 

Rf{n,t)«i, (B2) 

for any spatial angle Vt at any instant of time. The following argument is somewhat more 
intuitive. The deformation of the body is of the size i?/, and the external velocity changes at 
length-scales that are comparable with the correlation length ^. If the deformation is smaller 
than the correlation length (Fig. ^A) the external velocity does not change on the length- 
scale of the deformation, and the approximation is valid. On the other hand, if the correlation 
length is shorter than the deformation length-scale (Fig. |^B) the external velocities on the 
sphere and on the body are uncorrelated, and the approximation is unjustified. 
We turn to evaluate /. When i? << ^ it is clear that / can be made small by having 
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Vext small enough so that indeed Rf « C The more interesting case is R/^ > 1. The 
deformation / is determined by Eq. (D 

fl,m + Mfl,m + Ql,m = 0, (B3) 

where Qi^^ = ji[p ■ Vext\im- Clearly, 
An 

\Qlm\ < -^{\Ylni\)VexU (B4) 

where Vg^t is the typical magnitude of the external velocity. The average of the spherical 
harmonic is bound and of order one and therefore can be dropped off. Xifi^m is comparable 
with (eq. p3| ), therefore, in order that Rf « C, we must have Vext < ^iC,- A condi- 
tion that must be true for all values of / and especially for the smallest A; denoted A; rnin- 
Therefore, 

'^ext < A; min^ 

In most cases, however, we can find a stronger condition for the validity of the small de- 
formation approximation. We expect Qi^m to decline as the squared root of the number of 
independent surface elements, i.e. as ^/R, so that Xifi^m ~ ^^r- Therefore the condition, 
eq. (p^) , implies that Vext < XiR. Therefore, 

Vext < Xl rriinR- (B6) 

Both conditions can be easily maintained in a viscous fluid. The above conditions are 
general and depend on the specific system via A;. For example for a droplet with a surface 
tension energy and equal viscosities inside and outside, the minimal eigenvalue is A2 = 
1^ (where A is the surface tension constant and t] is the viscosity) and the condition is 
Vext < M^' '^^^^^ ^ viscosity much larger inside A/ = A2 = H and Vext < \^ 
where rj is the viscosity inside the droplet. 
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Figure captions 



Fig. 1 

The deformable body is described by a three dimensional scalar field The interior is 

the region where < 0, the exterior is the region where > and the outer surface of the 
body is the locus of the points obeying ip{r) = 0. 

Fig. 2 

The MSD for a fluid with a memory time scale. Fq = 
Fig. 3 

The Diffusion constant for typical separable random velocity correlations. Dq = 
Fig. 4 

The MSD for typical non-separable random velocity correlations with the scaling form, 
(t>{q-,t) — Cq~°'exp (^—Fqt^^, where we choose a = | and /3 = |. The MSD depends on the 
two non-dimensional variables: (^)'^^ and = CR°'~^'^'^r~'^ . The MSD scales for long 
times as in agreement with our scaling argument. 

Fig. 5 

The validity condition for the approximation. Fig. A: The deformation is negligible in re- 
spect to the velocity correlation length. Therefore the approximation holds. Fig. B: The 
deformation length is longer than the correlation length. The velocities on the surfaces of 
the sphere and droplet are uncorrelated. 
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FIG. 1. 
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FIG. 2. 
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